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Abstract 

Using the FOPI detector at GSI Darmstadt, excitation functions of collective 
flow components were measured for the Au+Au system, in the reaction plane 
and out of this plane, at seven incident energies ranging from 100^ MeV to 
800^ MeV. The threshold energies, corresponding to the onset of sideward- 
flow (balance energy) and squeeze-out eff'ect (transition energy) , are extracted 
from extrapolations of these excitation functions toward lower beam energies 
for charged products with Z >2. The transition energy is found to be larger 
than the balance energy. The impact parameter dependence of both bal- 
ance and transition energies, when extrapolated to central collisions, suggests 
comparable although slightly higher values than the threshold energy for the 
radial flow. The relevant parameter seems to be the energy deposited into the 
system in order to overcome the attractive nuclear forces. 

Keywords : Heavy ion collisions, nuclear matter expansion, sideward-flow, squeeze-out, 
radial flow, balance energy, transition energy. 

PACS numbers : 25.70.-z, 25.75.Ld 
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I. INTRODUCTION 



Collective motions of nuclear matter occurring in heavy ion collisions are of great interest 
since they are expected to provide information about the properties of hot and dense nu- 
clear matter and the underlying equation of state (EoS) [|I|. Flow effects were predicted by 
hydrodynamical calculations [|T],|2[ and experimentally evidenced at LBL-BEVALAC [0. At 
beam energies i?iab > 200A MeV, the interaction between nuclei is dominated by individual 
nucleon-nucleon scattering and the repulsive component of the mean field. This leads to a 
collective deflection of matter to positive angles in the reaction plane i.e. in the direction of 
the projectile remnants (sideward-fiow) . Conversely, at few tens of A MeV, the interaction 
is dominated by the attractive mean field, so that nucleons emitted in the reaction plane 
are deflected to negative angles 0]. At a certain intermediate incident energy, named the 
balance energy E-qal, the attractive component and the repulsive component of the interac- 
tions balance each other and consequently the flow crosses zero, changing from a negative 
sign at low energies to a positive sign at high energies. The balance effect was extensively 
investigated at GANIL 1^-0 and MSU [p|-pJ]l by measuring different colliding systems. It 



was also studied in the framework of theoretical models fT^-p^ . For semi-central collisions, 
the balance energy was found to be sensitive to the stiffness of the nuclear EoS and to the 
in-medium reduction of the nucleon-nucleon cross section \T3\. At lower impact parame- 
ters the balance energy is expected to be only sensitive to the in-medium nucleon-nucleon 
cross section ||T3| , pO| . Recent experimental results |p5|-p8l revealed also, in the same beam 
energy range, signatures of a change in the azimuthal emission pattern of mid-rapidity par- 
ticles, from an in-plane enhancement at low incident energies to the well known out-of-plane 
preferential emission, the so-called squeeze-out p9| , p0| , at higher energies. The incident en- 
ergy where this transition takes place (termed ii^TRA) is also found to be sensitive to the 
in-medium nucleon-nucleon cross section [p8| . 

The excitation functions of both sideward-flow and squeeze-out effects were measured 



with the FOPI detector |31| for the Au+Au system at seven incident energies between i^iab 
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100/lMeV and Eiab = SOOAMeV. The ability of the FOPI device to detect intermediate 
mass fragments {Z > 3) allows a cleaner identification of the collective flow signal and a 
better extrapolation of the measured excitation functions of the in-plane and of the out-of- 
plane flows toward their thresholds at low incident energies. We found however that the 
threshold energies are rather insensitive to the charge of the considered particles {Z > 1). 
We found also that the transition point seems to be located at higher beam energies than 
the balance point. On the other hand, the FOPI detector ensures a wide range of impact 
parameter collisions to be explored. This offers the possibility to investigate the centrality 
dependence of Ebal and -Etra, an aspect which was recognized to be crucial in this kind 
of study |[TT|JT^j23[| . This allows us to discuss for the first time this centrality dependence in 



conjunction with the threshold energy {E^iad) for the radial flow in highly central collisions. 
It shows that the three threshold energies {E^al, -^tra, and -Erad) niight be attributed to 
a common phenomenon, the relevant parameter being the energy deposited into the nuclear 
system. 

II. EXPERIMENTAL SETUP 

The data presented in this paper concern the Au+Au system at seven incident energies 
Ei.^^ = 100, 120, 150, 250, 400, 600 and SOOAMeV. They have been collected with the 



Phase I of the FOPI detector at the SIS/ESR accelerator facility, GSI Darmstadt. In 
its Phase I configuration, the FOPI detector covers in full azimuth the laboratory polar 
angles (Oiab) from 1.2° to 30°. It consists mainly of a highly segmented Forward Wall 
of plastic scintillators divided into two parts : the Inner Wall made of 252 trapezoidal 
scintillators which covers the Oiab range between 1.2° and 7.5°, and the Outer Wall made 
of 512 scintillator strips which covers the Biab domain from 7° to 30°. The Forward Wall 
provides an element identification and the velocity of the reaction products through energy 
loss and time-of-fiight measurements. Its segment structure allows to determine the velocity 
vector components. A complementary shell of 188 thin energy loss detectors is mounted 
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in front of the Forward Wall in order to achieve lower detection thresholds. This cluster 
detector is made of an ensemble of gas-filled ionisation chambers (Parabola) mounted in 
front of the Outer Wall, and thin plastic scintillator paddles (Rosace) combined with the 
Inner Wall. In order to reduce the background scattering in the air gas, a helium bag is 
placed between the target and the detectors. This setup measures simultaneously most of 
the light charged particles and intermediate mass fragments (up to Z = 15) emitted in 
the forward center-of-mass (cm.) hemisphere. Its high granularity allows high multiplicity 
events to be measured with a negligible multi-hit rate. The apparatus ensures a very good 
azimuthal symmetry which is an important feature for the study of the flow phenomena. 



III. EVENT CHARACTERISATION 

A. Impact parameter determination 

The measured events were sorted according to their degree of centrality using the stan- 
dard method based on the correlation between the multiplicity of emitted particles and the 
impact parameter. The multiphcity distribution of the charged particles detected in the 
Outer Wall exhibits the typical plateau for intermediate values, followed by a steep decrease 
at the highest multiplicities The highest multiplicity bin (named PM5) has been de- 
fined by cutting at half of the plateau value (the corresponding lower limits of the PM5 
multiplicity bin are given in Tab. 1). The remaining part of the multiplicity distribution 
has been subdivided into four equally spaced bins (named PMl to PM4) according to the 
procedure used by the Plastic Ball Collaboration . The results presented in what follows 



include only events belonging to the PM3-PM5 multiplicity classes where background con- 
tamination, estimated from measurements without target, is negligible. The mean impact 
parameter associated to each PM event class has been determined in the framework of the 
IQMD (Isospin Quantum Molecular Dynamics) model [|T6| , p4| by filtering the theoretical 



calculations with the realistic simulator of the FOPI detector. The resulting < b > values 
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and their r.m.s deviations are presented in Tab. || for the beam energies E'lab = 150 and 
400AMeV. Note that < b > values are the same within the r.m.s deviations at the other 
beam energies. It can be seen from Tab. y that the multiphcity criterium offers a large 



coverage of the impact parameter range. On the other hand, model studies |3^, using the 
so-called quality factor introduced by Cugnon and L'Hote allowed us to show that with 
the FOPI/Phase-I setup, the multiplicity criterium, as compared to other criteria, appears 
as the most appropriate one for exploring flow observables over a large impact parameter 
domain. 



B. Reaction plane reconstruction 

The reaction plane was reconstructed with the transverse momentum analysis devised 
by Danielewicz and Odyniec In order to remove autocorrelation effects, the azimuth of 
the reaction plane was estimated for each particle i in a given event as the plane containing 
the vector Qi and the beam axis, where Qi is calculated from the transverse momenta pi of 
all detected particles except the particle i : 

M is the multiplicity of the event and = 1 if ?/|°^ > 6,-1 if ?/|°'' < —6 and otherwise. 
y-^^ is the j*^ particle rapidity divided by the projectile rapidity in the cm. system. The 
parameter 6, choosen equal to 0.5, was introduced in order to remove mid-rapidity particles 
which have a negligible correlation with the reaction plane. According to [^], a boost 
velocity = p\/{nf^^ — m') {iv} is the mass of the particle i and Trf^^ is the sum of the 
projectile and target masses) was applied to each particle j in order to take into account 
the effects of momentum conservation due to the exclusion of the particle i. The influence 
of these effects on the observables considered in the present work was found to be of a few 
percent at the lowest beam energies (-Eiab = 100 and 120^ MeV) and negligible at higher 
incident energies. 
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The accuracy of the reaction plane reconstruction {i.e. the azimuthal deviation A$r 
of the reconstructed reaction plane with respect to the true one) was estimated for each 
event class by randomly dividing each event into two equal parts and by taking the one 
half of the angle between the Q vectors of the two subevents |^ . The azimuthal dispersion 



cr(A$R) was found to vary typically from ~ 20° to ~ 40° for the PM event classes under 
consideration. It is worth noting that the use of a heavy system such as Au+Au offers, 
in the beam energy range considered here, a good event characterisation both in centrality 
and reaction plane reconstruction as compared to lighter systems whose measured ejectile 
multiplicities are lower. 

IV. RESULTS AND DISCUSSION 
A. Sideward- flow 

The in-plane flow component (sideward-flow) is examined in terms of the normalized 
in-plane transverse momentum p^^ (p^^ = (Px/^) /Pcm where Px/A is the in-plane trans- 
verse momentum per nucleon and p^^^ is the projectile momentum per nucleon in the cm. 
system) as a function of the normalized cm. rapidity y^^^ {y^^^ is defined above). This 
normalisation, suggested in earlier works [^,|D|, is motivated by the fact that one obtains 



a scale invariant representation of the data in a fluid dynamical description of the collision. 
Figure |l] shows a typical example of p^^ (Fig. [l|.a) and < p^^ > (Fig. |l].b) versus y^^^ plots 
for Z = A particles detected in semi-central (PM4) reactions at i^iab = 250^4 MeV. Since the 
experimental apparatus covers only the forward cm. hemisphere, the plot has been mea- 
sured for positive y'-"^ rapidities and reflected for negative ones. As shown by Fig. |I].b, the 
dependence of < p^^ > on y^^^ exhibits the well known S-shape behavior ||3^ demonstrating 



the collective transfer of momentum between the backward and the forward hemispheres. 
The linear part of the curve in the participant region {i.e., at mid-rapidity) reflects the so- 
called side-splash effect while the fall-off starting just below the projectile rapidity (y^^^ = 1) 
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is caused by the bounce-off effect ET| . A quantitative measure of the amount of flow in the 



participant region of the collision is given by the so-called normalized flow parameter Fg°^ 
which is commonly defined as the slope of the < p^^ > versus y'^^'* curve at mid-rapidity 



d < > /d?/''°^ jjij ^. Technically the Fg parameter is obtained by fitting a 
polynomial function of the form : a + F^'^ x y^^^ -|- c x (y*-"^)^ to the data (Fig. |I[b). The fit 



was restricted to the linear branch of the S-shape curve. As shown in reference because 
the particle momenta are not projected onto the true reaction plane, their projections are on 
average biased downward by a factor 1/ < cos(A$r) > where A$r, as mentioned before, is 
the estimate of the azimuthal deviation of the reconstructed reaction plane with respect to 
the true one. The data shown in Fig. |l].b and all the Fg°^ values presented in what follows are 
corrected for this effect. The correction factors (1/ < cos(A$r) >) were typically ranging 
from 1.10 to 1.45 depending on the beam energy and the multiplicity bin. 

The precise evaluation of the acceptance effects on the in-plane flow is rather difficult due 
to the complexity of the different experimental constraints. This can only be investigated 
in the framework of realistic simulations where theoretical calculations are passed through 
the detector filter. In this context, we have used the IQMD model [0,0] which is known 
to reproduce quite well experimental flow data P5|-^. A few thousand of IQMD events 



were generated over a large range of impact parameters with the HM choice (Hard EoS 
plus a momentum dependent potential) of the nuclear interaction. This force is recognized 



as providing the best description of the observed trends in the in-plane flow data ||4^-^ 
For the present study, theoretical events were filtered applying geometrical cuts and en- 
ergy thresholds of the FOPI detector. They were presorted in accordance with the above 
mentioned procedure used for the data (see before). The in-plane flow was extracted with 
respect to the true reaction plane which is known in the model. Because of the limited 
statistics, apparatus effects could only be evaluated for light particles. We found that the 
experimental cuts only slightly affect < p^^^ > values in the forward cm. hemisphere. The 
observed deviations are mainly caused by the geometrical limit of detection at Giab = 30°. 
This cut biases down the F^^^ parameter by about 20% for Z = 1 and by less than 10% for 
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heavier particles in the PM5 event class at i?iab = 250 A MeV. Note that these effects decrease 



with increasing fragment size |^ because heavy particles, due to their low sensitivity to 



thermal fluctuations, occupy a smaller phase space. Therefore, in order to avoid misleading 
interpretations of the data, the F^^^ flow parameters presented in what follows include only 
the measurements of particles whose charge is > 2. On the other hand, it is worth noting 
that the effects of the Giab = 30° cut decreases with increasing impact parameters since 
peripheral and semi-central event topologies are less accentuated in the transverse direction 
as compared to central events. 

B. Squeeze-out 

The out-of-plane flow component was investigated from the azimuthal distributions 
dN/d^ ($ is the azimuthal angle of the detected particle relative to the azimuth of the 
reaction plane) around the beam axis of mid-rapidity particles, selected by imposing a ra- 
pidity cutoff -0.1 < < 0.1. It is now an established fact that the out-of-plane anisotropy 
increases strongly with the transverse momentum of charged particles p0|-|5^ . Therefore, in 
order to extract relevant information from the data, we determined the squeeze-out signal 
by choosing a pt cut which, within the acceptance, gives access to the largest momenta 
and sufficiently wide for statistics considerations. This pt window is 0.4 < p[^^ < 0.55, p[^^ 
being the particle transverse momentum per nucleon divided by the projectile momentum 
per nucleon in the cm. system. These rapidity and transverse momentum cuts used to 
extract the signal, define a portion of the phase space which is covered by the FOPI de- 
tector acceptance |51|. Figure |^ shows a typical dN/d^ distribution for Z = 3 particles 



in the PM4 event class at an incident energy of -Eiab = 250AMeV. A clear preferential 
emission is observed in the direction perpendicular to the reaction plane ($ = 90° and 
$ = 270°). This enhanced emission reflects the squeeze-out effect. The magnitude of the 
latter is commonly defined as the ratio i?N of the number of particles emitted perpendicu- 
lar to the reaction plane to the number of particles emitted in the reaction plane [00,0 : 
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= {N{90°) + N{270°))/{N{0°) + N{180°)). The i?N ratio is extracted by fitting a function 
of the form A^($) = + ai x cos($) + 02 x cos(2$) to the data (curve of Fig. g). Thus i?N 
is calculated as = (ao — 0.2) /{ao + 02). According to this definition, R^ < 1 and -Rn > 1 
are related to a preferential emission of matter in the reaction plane and out of this plane, 
respectively while R^ = 1 corresponds to a perfect azimuthally isotropic situation. The 
anisotropy ratio -Rn can be corrected, like the fiow parameter F^'^\ for the uncertainties due 
to fiuctuations of the reaction plane However, we found with the help of simulations 



that, for the low multiplicity events, the values of < cos^(A$r) >, which are the quantities 
involved in these corrections, were not determined with good accuracy. Thus, i?N ratios 
reported in what follows are not corrected for the effects of reaction plane fluctuations. Nev- 
ertheless, the possible influence of these effects on the observed trends will be discussed in 
the following. 

C. Balance energy 

Figure ^ shows the excitation functions of the scale invariant flow parameter F^^^ for 
different particles (Z = 2 to 5). A sudden decrease is observed in the incident energy 
region i?iab < 200AMeV. It is interesting to notice that this sudden change is much more 
pronounced in the case of the heavier fragments which are more sensitive to the collective 
motion. An extrapolation with Fermi functions allows us to estimate the balance energy 
(intersection with the abscissa) for different types of particles. Note that the balance energy 
values extracted from extrapolations of the data with other functions, such as logarithmic 



and second order polynomial ones, were found to be the same within error bars |^ . We have 
also verifled that the use of other scaling variables for the in-plane transverse momentum, 
such as < Px/pt > (as used in |Q) or < px > / < Pt >, leads to very similar results. 
The extrapolation for Z = 3 fragments leads to an intersection energy of -Ebal = 65 ± 
15AMeV for events corresponding to the PM4 bin. The resulting values for the other 
particles {Z = 2, A and 5) are the same within uncertainties as the one obtained for Z = 3 
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(see Tab. |T|). This confirms the observations estabhshed by studying hghter systems that 
the balance energy is independent of the size of the detected particle Our present 

balance energy point is somewhat larger than the one obtained for the same system from 



other experiments ||5^ , |55[] . This is probably due to the fact that our PM4 event class contains 
less central events than the one used in references Indeed, as it is shown in the 

following, the balance energy is found to decrease with decreasing impact parameters. With 
this in mind, the balance energy value extracted here is consistent with the systematics 
of the balance energy as a function of the mass of the combined projectile-target system 



obtained from MSU and GANIL results [Q. On the other hand, it is worth noting that for 
a heavy system like Au+Au, because of the strong Coulomb repulsion, the overall force is 
always repulsive p3[. Since the balance energy should correspond to the energy at which the 
attractive and repulsive component of the nuclear interaction balance each other, it must be 
evaluated without contamination of non nuclear contribution to the flow. For light systems 
this problem is less severe because of the much weaker Coulomb repulsion. Nevertheless, 
for heavy systems one may hope to extract the correct value of the balance energy by 
extrapolating the Fg°'* values from sufficiently high energies where Coulomb contribution is 
neghgible. 



D. Transition energy 

The dependence of the anisotropy ratio i?N on the collision impact parameter is presented 
in Fig. ^ at four incident energies going from 100 to 400^4 MeV. The signal includes here 
the contributions of all detected particles each weighted by its charge. By doing so we 
reconstruct a coalescence invariant quantity which makes meaningful the investigation of 
the anisotropy ratio as a function of the impact parameter and the beam energy. The 
geometrical impact parameter bg was obtained from the measured multiplicity distributions 
by assuming a sharp-cut-off approximation. This allows us to perform direct comparisons 
of data measured at different bombarding energies. As can be seen in Fig. ^, the correlation 
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between the i?N ratio and the impact parameter exhibits a very different trend as the incident 
energy decreases. At the highest bombarding energy (400AMeV) one can observe a bell- 
shaped distribution whose maximum is located at intermediate impact parameters (close to 
6/m). It is worth noting that the results obtained at higher energies are quite similar to 
those observed at i^iab = 400/1 MeV With decreasing beam energies, the shape of the 



correlation evolves gradually toward a different trend which is an evidence for a clear change 
in the emission pattern. Thus at -Eiab = lOOA MeV, with decreasing impact parameter, one 
observes a transition from a preferential in-plane emission (i?N < 1) to the squeeze-out effect 
characterized by an enhanced out-of-plane emission (i?N > !)• 

Before going to the interpretations of this behaviour, it must be pointed out that two 
effects might influence the ratio : i) the dispersion of the reconstructed reaction plane 
with respect to the true one, which tends to attenuate the magnitude of the signal and ii) 
the sideward-flow deflection which favors the emission of particles in the reaction plane. In 
both cases the magnitude of the effect is impact parameter and beam energy dependent. 
Therefore, in order to eliminate possible ambiguities in the interpretation of the experimen- 
tal observations in Fig. ^, it was necessary to examine the respective influences of these 
effects on the correlation between the anisotropy ratio and the impact parameter. A further 
complete analysis of the data has allowed us to show that the bell-like shape of the distribu- 
tion observed at high beam energies is preserved after taking into account both mentioned 
effects [^. On the other hand, the transition from _Rn > 1 to i?N < 1 at i?iab = lOOAMeV 
cannot be caused by one of these two effects. Indeed, fluctuations of the reaction plane 
tend to attenuate an anisotropy signal regardless of whether the ratio is larger or smaller 
than 1. Therefore, taking into account the corresponding corrections in Fig. ^, the transition 
effect would be even more pronounced. On the other hand, an extraction of the Rjss quantity 
around the flow axis would shift up the experimental points but at large 6g's, where the in- 



plane enhancement is observed, the flow angle is expected to be fairly low I^S] in particular 
at the lowest beam energy (-Eiab = lOOAMeV) which is close to the balance energy (Fig. 
Let us now go back to the interpretations of the experimental observations of Fig. 
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At high energies, the maximum located near bg = 6/m is consistent with an expansion- 
shadowing picture, i.e. an expansion of the compressed matter in the central region of the 
collision which is hindered by the presence of cold spectator remnants. It is worth noting 
that recent IQMD calculations for neutrons |Q predict a similar bell-shaped correlation 
with a maximum around 7/m. At low incident energies (-Eiab < ISOAMeV), a clear evidence 
for a transition from an enhanced in-plane emission pattern to a preferential out-of-plane 
emission is observed when bg decreases. The results show that this transition takes place 
close to Eiah = lOOA MeV for collisions with impact parameters bg ~ 6/m. This effect was 
already observed for a lighter system pSf and very recently for the same system Au-|-Au p8 



It might be attributed to a change from a collective rotational behaviour governed by the 
attractive mean field at low energies, to the high energy squeeze-out effect resulting from 



the repulsive pressure built up during the high density stage of the collision [25,28 



The transition energy -Etra, corresponding to an azimuthally symmetric distribution 
(i?N = 1), can be evaluated from the excitation function of the anisotropy ratio -Rn- This 
has been already investigated in a previous analysis where the experimental -Etra values 
were compared to the predictions of the IQMD model. The excitation functions of the i?N 
ratio are presented in Fig. ^ for Z = 2 and 3 particles emitted in semi-central (PM4) 
collisions. Since events were selected here over a large multiplicity bin, the large statistics 
allowed us to investigate the anisotropy signal in a reduced high pt window as compared 
to the previous one. This pt condition (fixed as 0.5 < p[^^ < 0.55) was chosen in order to 
extract from the data the largest -Rn magnitudes within the acceptance of the detector. As 
one can see from Fig. ^, the ratio tends to saturate above i^iab = 250^4 MeV and seems 
possibly to decrease at higher energies. On the other hand, the behaviour of the signal at 
low beam energies exhibits the same sudden change as in the case of the excitation functions 
of the sideward flow. 

In order to extract the transition energy, the data points in the beam energy range 
100-400^4 MeV were fitted with a Fermi function (curves of Fig. ^. The values of -EtrAj 
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reported in Tab. were determined for Z = 2 and 3 particles at the intercept of i?N = 1 
with the curves. Note that the intersection energies were found to be quite insensitive to 
the form of the fitting function. They seem to be, as observed in the case of the balance 
energy, independent of the type of the detected particle within the error bars. It can be seen 
from the results presented in Tab. ^ and Tab. ^ that the transition energy for the PM4 
event class (-Etra — lOOAMeV) is somewhat larger than the corresponding balance energy 
(E'bal — 65AMeV). This finding agrees with recent theoretical calculations for the Ca+Ca 



system ||2^. It could be explained by angular momentum effects or by the fact that at 
the balance point the compression is not high enough to generate an enhanced out-of-plane 
emission. 

On the other hand, the i?N ratio, which increases with the particle transverse mo- 
mentum for i?iab > ISOAMeV [^0|-|52|, has been extracted over a narrow pt window 
(0.5 < < 0.55). This could suggest that the extracted transition energy is pt de- 
pendent. Nevertheless, as shown by the insert in Fig. ^, at E\^\^ = lOOAMeV where the 
transition effect takes place, the ratio remains nearly close to 1 in the whole explored 
p[^^ domain. Finally, as mentioned before, the i?N ratios are not corrected for fluctuations 
of the reconstructed reaction plane because the < cos^(A$r) > values cannot be accurately 
determined. Therefore, in order to estimate how the transition energy could be influenced 
by these fluctuations, we have extracted -Etra from the values obtained by correcting the 
measured anisotropy ratios using the < cos^(A$r) > factors calculated for flltered IQMD 
theoretical events, where the true reaction plane is known. The resulting transition energy 
value, for Z = 3 particles, was found to be -Etra = 111 ± lOAMeV in the PM4 event class. 
This value is very close to the one obtained without correction (Tab. |rv|) . 

E. Centrality dependence of threshold energies 

The balance and transition energies extracted from extrapolations of the measured Fg^-* 
and i?N excitation functions in the three different PM multiplicity bins are shown in Fig. as 
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a function of the collision centrality. Ebal and -Etra are expressed here not in terms of the 
projectile energy but of the corresponding energy in the cm. system. They were obtained, 
due to considerations of detector acceptance effects and statistics, from the measured exci- 
tation functions of Z = 3 for -Ebal and Z = 2 for -Etra- As it can be seen, the transition 
energy is larger than the balance energy over the whole explored impact parameter range. 
Both threshold energies -Ebal and -Etra increase with increasing impact parameters. Such a 
behaviour has been also very recently observed in the case of the balance energy of the Ar+Sc 
system |]ll| and is predicted by the IQMD model for the Ca+Ca system [^. This increase 



of -Ebal and -Etra with b may indicate that the threshold of flow effects is related to the en- 
ergy deposited locally into the overlap zone of the collision : if the local temperature is large 
enough, the generated pressure can overcome the attractive nuclear forces. When going from 
peripheral to central collisions, the deposited energy becomes larger and the onset of flow is 
therefore expected to take place at a lower incident energy. At bombarding energies close to 
the Fermi energy, the participant picture is not well developed like at higher energies and 
the local temperature depends sensitively on the local heat relaxation-time and on collision 
time. With this in mind, it is interesting to compare the threshold energies in finite impact 
parameter collisions with the threshold energy for the radial expansion in central collisions. 
The latter was determined in reference |^9l| by extrapolating the excitation function of the 



mean radial flow velocity, measured in central collisions, toward low incident energies. Its 



value was found to be about -Erad = 35 ± lOAMeV |59| which correponds to a cm. energy 
of 8.7 ± 2.5AMeV. Note that this value is consistent with recent results for the Au+Au 
reaction obtained at -Eiab = 35AMeV [Q. Now, considering again Fig. ||, a rough linear 
extrapolation of our balance and transition energy points toward < b > = leads to values 
of 12 ± 6y4MeV and 16 ± 5^4 MeV, respectively. These values are close to each other and 
consistent with the value of -Erad within the error bars. It is still premature to go to more 
quantitative interpretations because of the large experimental uncertainties. Nevertheless, it 
is tempting to speculate on the basis of the observations in Fig. ^ on a common phenomenon 
at the origin of radial flow, sideward-fiow and squeeze-out effects which takes place for the 
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same deposited energy : when nuclear matter is heated above a certain hmit, the attractive 
nuclear forces are counterbalanced by the thermal pressure and the system starts to expand 
above this threshold at which the overall force becomes repulsive. This expansion of nuclear 
matter tends to an azimuthally symmetric pattern (radial flow) in highly central collisions, 
while in semi-central collisions the expanding matter is pushed to the side (sideward-flow) 
and the presence of cold spectator remnants hinders participant nucleons to escape in the 
reaction plane, which gives rise to a preferential out-of-plane emission (squeeze-out). This 
scenario for the squeeze-out effect is consistent with the observed correlation between the 
anisotropy ratio and the impact parameter (Fig. Differences on the various threshold 
energies can be expected on the bases of effects like collision geometry, system size. Coulomb 
contribution, angular momentum dissipation, non equilibrium dynamics, etc. On the other 
hand, it is worth noting that the cm. energy of 8.7AMeV where the radial flow sets in, is 
close to the one of the caloric curve, reported by the ALADIN Collaboration ||61|j6^ , where 
the temperature grows strongly again after the presently debated plateau. As pointed out 
in ref. |HD|, this observation suggests that the onset of flow could be interpreted as a possible 
signature of a liquid-gas phase transition. Indeed, as discussed before, the onset of flow 
phenomena indicates a change in the reaction scenario from a global repulsive mechanism 
to a global attractive mechanism when decreasing the incident energy. In a fluid dynamical 
vision of the collision, it could also be seen as a manifestation of a liquid-gas phase transition 
since this leads to a sudden decrease of the repulsive pressure and, consequently, of the flow 
magnitude. 



V. CONCLUSION 

Collective effects of nuclear matter in Au+Au collisions at incident energies ranging 
from lOOA MeV to SOOA MeV were measured for light particles and intermediate mass frag- 
ments with the FOPI detector at GSI. The centrality and beam energy dependence of both 
sideward-flow and squeeze-out effects were investigated. The evolution of the squeeze-out 
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magnitude i?N with the impact parameter is found to change drastically with the incident en- 
ergy. At high incident energies the squeeze-out signal exhibits a bell-shape with a maximum 
located at intermediate impact parameters. This trend is consistent with an expansion- 
shadowing scenario where the expansion of highly compressed participant nuclear matter 
is hindered by the presence of cold spectator remnants. At low energies, the data show 
clearly a transition from an in-plane preferential emission to an out-of-plane enhancement 
when the centrality increases. This phenomenon might be attributed to a change from a 
collective rotational motion for large impact parameters to the squeeze-out effect for smaller 



impact parameters [^,^. The balance and transition energies, corresponding to the onset 
of sideward-fiow and squeeze-out effects, respectively, have been evaluated from extrapo- 
lations toward lower beam energies of the excitation functions of the scale invariant flow 
parameter F^^^ and the anisotropy ratio R^. Both of them are found to be, within errors, 
independent of the size of the detected particle. The transition energy is larger than the 
balance energy. The extrapolation of the centrality dependence of both threshold energies 
toward 6 = leads to a value which is close to the threshold energy for the radial expansion 
in central collisions. This suggests that the same phenomenon could be at the origin of the 
three processes. The relevant parameter seems to be the energy deposited into the system in 
order to counterbalance the attractive nuclear forces although compression effects can also 
be present. In central collisions the repulsive pressure is expected to be the highest and the 
participant matter can expand freely in all directions, while for 6 7^ 0, besides complex geo- 
metrical and dynamical effects, the presence of the spectator matter causes the appearence 
of sideward-fiow and squeeze-out. Comparisons of the present experimental results with 
transport model predictions should provide interesting information about the in-medium 



reduction of the nucleon-nucleon cross section 13,20,28 
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TABLES 



TABLE I. Lower limit PM5 y of the PM5 multiplicity bin at the different beam energies. 



Eiab (^MeV) 


100 


120 


150 


250 


400 


600 


800 


PMSj 


28 


31 


36 


44 


55 


62 


70 



TABLE II. Mean impact parameter < b > for each experimental PM multiplicity bin at incident 
energies of 150 and 400A MeV. < 6 > is determined in the framework of the IQMD model with a 
hard momentum dependent interaction. Errors represent the r.m.s. deviation of the b distributions. 



Multiplicity bin 


PM3 


PM4 


PM5 


Eiab=150^MeV 


8.8 ±1.9 


5.0 ±1.9 


3.4 ± 1.4 




7.0 ±1.1 


4.1 ±1.5 


3.2 ± 1.1 



TABLE III. Balance energy for Z = 2 to 5 particles under the PM4 multiplicity cut (see text). 
Errors correspond to systematic uncertainties. 



z 


2 


3 


4 


5 


Ebal (A MeV) 


56.0 ±21.4 


65.2 ± 14.9 


68.3 ±11.1 


64.5 ± 15.1 
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TABLE IV. Transition energy for Z = 2 and 3 particles under the PM4 multiplicity cut (see 
text). Errors correspond to statistical uncertainties only. 



z 


2 


3 


£;tra (AMeV) 


98.9 ± 7.6 


107.0 ± 6.9 
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FIGURES 

FIG. 1. Upper panel : Normalized in-plane transverse momentum versus the normal- 
ized cm. rapidity for Z = 4 particles detected in semi-central (PM4) Au-|-Au collisions at 
Eiah = 250AMeV. The plot is obtained by assuming a forward/backward symmetry. The different 
grey levels correspond to different linear cuts in multiplicity. The 9c.m. = 30° cut is represented 
by the solid white curves. 

Lower panel : Mean normalized in-plane transverse momentum < > as a function of y^^^ for 
Z = 4 particles detected in semi-central (PM4) Au (250^4 MeV)-|-Au colhsions. The data (open 
stars) are measured only for y^'^^ > and then reflected around the origin (full stars). Data points 
are larger than the corresponding statistical uncertainties. The solid curve is the result of the fit 
described in the text. 

FIG. 2. Azimuthal distribution of mid-rapidity (—0.1 < y(°) < 0.1) Z = S particles measured 
in semi-central PM4 collisions at -Eiab = 250A MeV. $ is the particle azimuthal angle around the 
beam axis, with respect to the reaction plane. The distribution is extracted with the transverse 
momentum cut 0.4 < p[^^ < 0.55. The solid curve is the result of the fit described in the text. 
Error bars correspond to statistical uncertainties. 

FIG. 3. Excitation functions of the normalized flow parameter F^^^ measured for different 
particles [Z = 2 to 5) in semi-central PM4 collisions. The values are corrected for fluctuations 
of the estimated reaction plane. The solid lines correspond to the fits to the data points from 
E\ah = 100 to 400^4 MeV with Fermi functions. Error bars correspond to systematic uncertainties, 
estimated to 20, 20, and 10% for Ei^^ = 100, 120, and 150^ MeV, respectively and less than 10% 
for E'lab > 250A MeV. These errors are larger than statistical uncertainties. 

FIG. 4. Anisotropy ratio i?N as a function of the geometrical impact parameter bg for Au-|-Au 
collisions at Ei^h =100^ MeV (crosses), 150AMeV (triangles), 250AMeV (squares) and 400AMeV 
(circles). i?N includes the contributions of all detected particles with 0.4 < pf'^ < 0.55, each being 
weighted by its charge. Error bars correspond to statistical uncertainties. 



25 



FIG. 5. Excitation functions of the ratio for Z = 2 (crosses) and Z = 3 particles (triangles) 
detected in semi-central PM4 collisions under the condition 0.5 < p^^^ < 0.55. The curves represent 
the results of a fit with a Fermi function to the data points from £Jiab = 100 to 400^ MeV. The 
insert shows the pf'^ dependence of the i?N ratio at an incident energy of lOOA MeV. Error bars 
correspond to statistical uncertainties. 

FIG. 6. Impact parameter dependence of the cm. threshold energies -Eral (circles) and -Etra 
(triangles), determined from the measurements of Z = 3 and Z = 2 particles, respectively. < b > 
has been determined within simulations using the IQMD model at Eia,h = 150^4 MeV (see Tab. 2). 
The dotted lines are linear fits to the data. Horizontal error bars represent the r.m.s. deviation 
of b distributions. Vertical error bars correspond to systematic and statistical errors for -Ebal and 
Etra, respectively. 
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